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We consider synchronization of coupled chaotic systems and propose an adaptive strategy that 

(N ' 

— ■ aims at evolving the strength of the coupling to achieve stability of the synchronized evolution. We 

<D ■ 

. test this idea in a simple configuration in which two chaotic systems are unidirectionally coupled 

(a sender and a receiver) and we study conditions for the receiver to adaptively synchronize with 
the sender. Numerical simulations show that, under certain conditions, our strategy is successful 

'j-l', 

D , in dynamically evolving the coupling strength until it converges to a value that is compatible with 

o 



-a 

O 

o 



synchronization. 



I. INTRODUCTION 

It is known that stability of the synchronized solution of coupled dynamical systems depends on the strength of the 
interaction. Motivated by the observation that real-world networks are characterized by evolving adapting connections, 
some studies have started to take in consideration how the time variability of the strength of the couplings may affect 



> 

the synchronization dynamics. A few papers have appeared on adaptive strategies, in which the coupling strengths 

(N ; 

, are dynamically adjusted based on information on the states of the systems, e.g., in order to enhance or guarantee 

o; . 

synchronization. In this paper we address the problem of how to evolve the coupling strength from a given initial 

o ■ 

condition to make it converge to a value which corresponds to stability of the synchronized evolution. 



Previous works have addressed adaptive synchronization of chaos. For instance, in [l|, Q], a problem was studied 



in which adaptation was needed in order to achieve synchronization in the presence of external unpredictable events 
affecting the communication between the coupled systems. As an applicative example, consider the case of two sensors 
that seek to synchronize through a continuous signal that they exchange via wireless communication; if an object 
moves across the communication pathway between the sensors, this causes an attenuation on the strength of the 
received signal, which may disrupt synchronization. Therefore, in [l, 2|, adaptive strategies were introduced, having 
the aim of maintaining synchronization with respect to such external (unpredictable) perturbations. 

Here, we will focus on a different but related problem, namely we consider a simple configuration in which a sender 



2 

is connected to a receiver and communication between the two is such that attenuations affecting the received signal 
are negligible; the problem is rather for the receiver to adaptively choose an appropriate coupling form (in what follows 
we will specify our proposed problem to be that of choosing an appropriate coupling strength) with the received signal 
in order to achieve synchronization with the sender. 

To better illustrate this problem, we introduce the equations for the sender and the receiver, 

ii(t) = F(x!(t)), (la) 
x 2 (t) = F(x 2 (t)) + j(t)(H( Xl (t)) - H(x 2 (t))), (lb) 



where x\(t) (^(i)) is the m-dimensional state of the sender (receiver) system, Xi = [xn, Xi2, Xj m ]; F(x) is the 
dynamics of an uncoupled system (hereafter assumed chaotic), F : R m — > R m ; H{x) is an output function, H : R m — > 
R m ; and j(t) is a time varying scalar function measuring the strength of the coupling. Here we assume that H{x) can 
be rewritten as H(x) — Hh(x), where H — [Hi,H.2, ■■■,T~t m ] is a constant m-vector and h : R m — > R is a scalar output 
function. In many practical situations, communication between the connected systems involves only a subset of the 
dynamical state variables of the systems. Therefore, in this paper we consider a situation in which synchronization 
has to be achieved based solely on the scalar signal, h(xi(t)), that system 2 receives from system 1. Our proposed 
problem is to devise an adaptive strategy to dynamically evolve the coupling strength j(t) from an arbitrary initial 
condition 7(0) to a value that is compatible with synchronization. 

As an example of an application for our proposed strategy (to be specified in Sec. II), we consider a sender-repeater 
communication scenario, in which the first device, say a wifi antenna, sends a signal, and the second identical device 
needs to reproduce precisely the same signal to other systems that cannot get the original one. As another example, our 
strategy could be used in applications that take advantage of synchronization of chaos for identification and prediction 
of the dynamics of unknown real systems 3|, U, |5|, |6|. In particular, for the problems addressed in 3|, |4|, |5|, |6j, the 
receiver is a model system which is meant to replicate the dynamics of an unknown true system (the sender) and a 
technique is proposed that aims at identifying the true system parameters by making the receiver synchronize with 
the sender. The technique relies on the choice of a coupling strength which is compatible with synchronization, but for 
this case, an extra difficulty is due to the fact that the dynamical form of the sender's F(x) is unknown at the receiver 
and therefore one is unable to compute the 7-range of stability (i.e., compute the master stability function; for more 
details on this subject the reader is referred to the background section) beforehand. Therefore, if one has to deal with 



such a problem and unless one wants to exclusively rely on a try-and-error approach, it becomes necessary to introduce 
an adaptive strategy to dynamically evolve the coupling strength. We anticipate that the adaptive strategy that we 
present here does not use information from the master stability function and is based exclusively on minimization of 
the squared synchronization error [h(xi(t)) — h(x2(t))] 2 . Therefore, we envision our strategy would be particularly 
useful in applications that use synchronization as a tool to achieve parameter estimation and prediction of unknown 
systems (e.g., for all the applications described in [3j, \4, la, la] ) • 



A. Background 



We note that Eqs. (|la|) and (|lbj) admit a synchronous solution of the type, 



x 2 {t) = Xx{t), 



(2) 



which obeys the sender evolution (JTaJ) . Our proposed problem is to devise a strategy to evolve "f(t) from a given 
initial condition 7(0) in order for the receiver system to synchronize with the sender. Let us now assume that "/(t) 



is constant and equal to 7. Following 7], la llA llll , we know that depending on the choice of the functions F and H, 
synchronization is possible in a certain range of values of 7, say T . To see this, let us linearize the receiver system 
equation 11b[) about @, obtaining, 



5± 2 {t) = [DF(x 1 (t))-*yDH(x 1 (t))]6x 2 (t), 



(3) 



where DF(x\{t)) and DH(x\{t)) are the Jacobians of the functions F and H evaluated about the synchronous solution 
([2]) . Then stability of the synchronous solution depends on the maximum Lyapunov exponent associated with Eq. ([3]) 



(this is a classic result, see 



m 



10( for more details). Once the functions F and H are given, it becomes possible 
to study the dependence of the maximum Lyapunov exponent of (j3]) on 7. The function that associates 7 with the 
maximum Lyapunov exponent of ^ is usually called a master stability function 9|. 



From 



Kj ] and subsequent papers, we know that the range of values of 7, say T, in which the master stability 
function associated with Eq. ([3]) is negative, can be either bounded or unbounded. An example of a bounded range for 
7 is reported in Fig. 1, where for F(x) being the equation of the Rossler system, m = 3, and hixi) = xn, TL = [1, 0, 0], 



the master stability function is found to be negative in the range T = [0.3, 5]. Moreover, we observe from Fig. 1 that 



the master stability function has a minimum at 7 ~ 1.7 (corresponding to the maximum rate of contraction for ([3]) 
and therefore the fastest possible convergence to the synchronized evolution) . Another example of a bounded stability 
range T is reported in [10( for F(x) being the equation of the Lorenz system, m = 3, and h{xi) — x&, TL = [0, 0, 1]. 
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Y 



FIG. 1: Master stability function corresponding to Eq. (3), for our choices of F(x) being the Rossler equation and 
h(xi) — xn, TL — [1, 0, 0]. The zero ordinate line is shown as a guide to the eye (dotted line). 



The problem we address in this paper (a strategy to evolve the s 



■trength of the coupling for 



synchronization) has already been addressed in a number of papers 3 
following adaptive strategies were proposed, 



11 



12 



Q 



achieving stable 



Namely, in 



11 



12, lia], the 



7(t) = a\h( Xl {t)) - h(x 2 {t))\, 
\h( Xl (t)) -h(x 2 (t))\ 



7(*) = 



l + \h(x 1 (t))-h(x 2 (t))\' 



in [12,13], 
in [11], 



(4a) 
(4b) 



where a > is a suitable scalar gain. Equation (|4b[) is approximately equal to Eq. (|4a| for \h(xi(t)) — h{x 2 {t))\ <C 1, 
and is approximately equal to j(t) — a, for \h(x±(t)) — h(x2(t))\ ^> 1. It is worth noting that both strategies are 
based on increasing the coupling strength 7 proportionally with the absolute value of the synchronization error (with 
a saturation in the case of l[4b])). 



We note that the strategies proposed in 



11 



12L 1X31 ] may not be the best suited to deal with coupled dynamical 
systems whose stability is governed by a master stability function as that shown in Fig. 1. For example, for the case 
shown in Fig. 1, if 7(0) is larger than 5, then increasing 7 will not be effective in order to achieve synchronization. 
Moreover, even if 7(0) € T, strategies (|4al4bj) may eventually increase 7 outside the range T (e.g., beyond 5 in Fig. 
1) and yield desynchronization even in a case in which the systems would have synchronized with no adaptation (i.e., 
with 7 = 0). 
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Another adaptive strategy has been proposed in 3[, based on the following adaptation for 7(f), 

7 (t) = -07(f) + - M«a(*))) 2 ]. (5) 

where a > and 5(2;) is a function such that g(x) ~ x for small x and is bounded by a constant C for large x. We 
note from ((SJ, that for large t, 7(f) converges to a _1 <7[(ft.(a;i(f)) — /i(x2(f))) 2 ] and therefore for a monotonely increasing 
g(x), 7(f) increases monotonely with the squared synchronization error. Hence, ([5]) is subject to the same sort of 
limitations as (|4a|4b|) when stability of the synchronized evolution is described by a master stability function as that 
shown in Fig. 1. 

In what follows, we will present a novel adaptive strategy, which is aimed at evolving 7(f) from an arbitrary initial 
condition 7(0) to converge on T. We note that the following limitations may apply to some of the previously reported 
strategies (|4a|4b|5j) . 

1. The strategies (|4al4bl5l) are based on increasing the coupling strength with the synchronization error. Yet, for 
the case of a bounded range of stability T, either increasing or decreasing the coupling strength 7 may be needed 
in order to achieve synchronization. 

2. The strategies (|4al4bl5[) are dependent on the choice of the initial condition 7(0). 

3. The strategies (|4al4bl5[) do not seek to make the coupling 7(f) converge in a neighborhood of its optimal value 
(e.g., the minimum of the master stability function which is observed to be at about 7 = 1.7 in Fig. [T]). 

In what follows, we focus on the common situation that the 7-range of stability is bounded from both below and 
above (e.g., as shown for the case of the Rossler systems in Fig. [TJ. We devise an adaptive strategy that is based 
on both increasing/decreasing 7(f) in order to reach synchronization and will be shown to be independent of the 
initial condition 7(0), as far as 7(0) is not too distant from T. Moreover, though our problem is characterized by a 
continuous of solutions (represented by the range T), our hope is that 7(f) will converge not too far away from the 
value of 7 = 1.7 which corresponds to the minimum of the master stability function, that is to the maximum rate of 
contraction towards the synchronization manifold. 

We wish to emphasize that our strategy described here is only a possible alternative approach to others already 



presented in the literature (e.g., [3 



11 



12 



131 ]) and we assume that many other possible solutions can be found to 
the general problem addressed in this paper. At the same time, we hope our attempt will motivate further studies in 
the same direction, which will provide better and better solutions to our proposed problem. 
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The rest of the paper is organized as follows. In Sec. II, we present our adaptive strategy. In Sec. Ill, we show the 
results of numerical simulations involving our adaptive strategy. In Sec. IV, the conclusions are presented. 



II. ADAPTIVE STRATEGY 

In this section we present our adaptive strategy. We introduce a potential/cost function ^(t), 

*(i)=< [h{ Xl (t)) - h(x 2 (t))} 2 >„, (6) 
where < G(t) > v denotes the sliding exponential average 

j* e -Ht-t')G(t')dt'. From ©, we note that #(t) > 0; 
moreover, due to the chaotic nature of the x's, ^(t) =0 can only be realized if x\{t) and X2{t) are in the synchronization 
manifold, i.e., X\(t) = X2(t). Thus we seek to evolve ^(t) in order to minimize ^(t). To this aim we introduce the 
following gradient descent relation, 

tf . _„« _ „ < lh(xi(m _ >„= WG. (7) 

where (3 > 0, and 

^ = -vG{t) + [h( Xl (t)) - h{x 2 {t))]^^- (8) 

We are interested in how evolves in time. To this aim, we write, 

d dx2 dx2 r dF dH.dx2 rTT/ , xx . , ... , . 

To conclude, our adaptive strategy is fully described by the following set of differential equations, 



7 = 20G, (10a) 
dh T 

G = -vG(t) + [h(x 1 (t))-h(x 2 (t))]— £, (10b) 

0x2 

* = te _ W e + [jff(:Ei W) _ ff(a;2( * ))] - (10c) 
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III. NUMERICAL EXPERIMENTS 

In this section we present numerical experiments involving our proposed adaptive strategy (|10p . We integrate the 
set of Eqs. (fTj) . ([TO]) , and we specify our sender-receiver dynamical systems to be described by the Rossler equation, 
in = 3, 



F(xi) = 



— X i2 — Xtf 

xn + 0.2x i2 
0.2 + (x a - I0)x a 



(11) 



1, 2, coupled through h(xi) — xn, TL — [1, 0, 0]. Fig. [2] shows the average synchronization error E defined as, 



E 



(*a -h)pj tl 



\h(x 2 (t)) - h(xi(t))\dt, 



(12) 



from integration of Eqs. Jl} with j(t) = 7, versus 7; the normalization factor p = < (h(xi)— < h(xi) >) 2 > x / 2 , with 
< ... > indicating the time average, is calculated from dynamics of system in the synchronous state, i.e., using 
dynamics from Eq. (|la|) . 




FIG. 2: The average synchronization error E from integration of Eqs. |T} versus (constant) 7, for 7 in the range [—1, 10], 
ti = 0.9 x 10 2 , t 2 = 10 2 , pi = 7.45. 



We first consider the case that 7(0) is not far away from the stability range T = [0.3,5]; namely we consider two 
cases: in the first case we choose 7(0) = 0, and in the second 7(0) = 5.3. In our numerical simulations, system (fl~a|) 
and system (| 1F>|) are evolved from random initial conditions on the Rossler attractor. We choose v — 0.1 so that 
is larger than the characteristic time scale of a chaotic oscillation T ~ 2ir, and (3 to be small enough that 7 changes 
slowly on the time scale of the chaos, that is {3 = 10~ 4 . 
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The results are shown in Fig. [3] (for the case 7(0) = 0) and Fig. [4] (for the case 7(0) = 5.3). In both cases, we see 
that for large enough t, 7(f) converges to T = [0.3,5], e.g., it approaches 2.32 for the case in Fig. [3] and 1.78 for the 
case in Fig. [4] These values are not far away from the minimum of the master stability function at 7 ~ 1.7 (see Fig. 
1). 




_ 2 o I 1 1 1 1 1 

O 20 40 60 80 100 




O 20 40 60 80 100 

t 



FIG. 3: The upper plot shows the time evolution of Xn(t) (dashed line) compared to X2i(t) (continuous line), while the lower 
plot shows the time evolution of j(t). 7(0) = 0, (3 — 10 -4 , v — 10 -1 . 
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FIG. 4: The upper plot shows the time evolution of xu(t) (dashed line) compared to xn(€) (continuous line), while the lower 
plot shows the time evolution of 7(f). 7(0) = 5.3, /3 = 1CP 4 , v = 10 _1 . 

When 7(0) is distant from the synchronization range T, our proposed problem becomes more complicated. In fact, 
since it may take a long time for j(t) to reach T, X2(t) may move away from x\{t) and the master stability function 
analysis, which describes local stability about the synchronization manifold, may not apply anymore. 

Figure [5] shows a case in which j(t) is evolved from 7(0) = 6.5, for (3 = 1CP 4 and v = 1CP 1 . For this case, 
synchronization is not achieved, and the emergence of an unexpected phenomenon is observed (namely, the emergence 
of a new attractor, corresponding to X2i{t) ~ 0, with both |x23 | and |x22(£)| growing in time in such a way that 
X23(t) — —X22(t)). In particular, from Fig. Owe see that j(t) first decreases from 7(0) = 6.5 to a minimum value of 
about —1.5 and then increases again up to about 44; it is worth noting that 7(i) crosses twice the synchronization 
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range T = [0.3, 5] and neither the first nor the second time it converges to it. We have repeated the experiment in 
Fig. [5] several times, eventually observing a different behavior, that is, j(t) increases indefinitely above 7(0) = 6.5 
and never goes through the synchronization range T (not shown). 

Based on the above observations, we conclude that, in order to improve the effectiveness of our strategy, we need 
to appropriately tune the rate of change of "f(t). In particular, we note that (i) in the case that j(t) is close to T, 
our strategy would benefit from bounding jit) so as to allow enough time for "/(t) to converge in T; (ii) in the case 
that "f(t) is distant from T, our strategy would benefit from choosing large values of (3 so as to quickly adjust ^(t) to 
converge towards V . 
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FIG. 5: The upper plot shows the time evolution of xn(t) (dashed line) compared to X2i(t) (continuous line), while the lower 
plot shows the time evolution of 7(i); in the lower plot, the zero ordinate line is shown as a guide to the eye (dotted line). 
7(0) = 6.5, P = 1CT 4 , v = HT 1 . 



Therefore, we propose to replace -y(t) — 2(3G(t) in (|10a[) by, 



where Eq. (fT"3|) is approximately equal to (|10a[) in the case that \G(t)\ <C 1, while is approximately equal to -y(t) = 
2/3 x sgn(G(t)) in the case that \G(t)\ 3> 1. Note that Eq. (1 13[) corresponds to adding a saturation on 7(f), which 
now is constrained to approximately lye in the range [—13, 0\. 

Figure [S] (Figure [7]) shows the results of numerical simulations in which we have tested our modified adaptive 
strategy, given by Eqs. (fTj). (fT3)) . (|l0b[) , (|10c|) , for a case in which 7(0) = —1 (7(0) = 10). For our experiments in Figs. 
|5]and [71 we have chosen j3 = 10 _1 . As can be seen, in both cases 7(t) is observed to converge to the synchronization 
range T, though it starts from initial conditions that are far away from it. 
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FIG. 6: Modified adaptive strategy. The upper plot shows the time evolution of xu(t) (dashed line) compared to X2i{t) 
(continuous line), while the lower plot shows the time evolution of y(t). 7(0) = — 1, /3 = lO -1 , v = 1. 
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FIG. 7: Modified adaptive strategy. The upper plot shows the time evolution of xii(i) (dashed line) compared to X2i{t) 
(continuous line), while the lower plot shows the time evolution of j(t). 7(0) = 10, (3 = 10 -1 , v = 1. 

A. An experiment with the Lorenz system 

We now consider an example in which our sender-receiver dynamical systems arc described by the Lorenz equation, 
to = 3, 



l0(Xi2 - Xn) 

23x l i - x i2 - XnXiz 



(14) 



n 



i = 1, 2. From [14| we know that when two Lorenz systems (|14l) are coupled as in Jl} with constant coupling j(t) = 7, 
h(xi) = Xi3, TL = [0,0,1], synchronization is stable in a bounded interval of the coupling strength 7, that is for 

7 e r = [1.2,6.5] . 

Figure [8] shows the average synchronization error E, defined in Eq. (|12p. for two Lorenz systems coupled as in 
Eqs. fl}, with constant coupling 7(f) = 7, h(xi) — 2^3, TL = [0,0, 1], as function of the parameter 7; for this case, we 
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17], which 



observe that the transition to the non-synchronous state is characterized by on-off intermittency [15j, 116 
makes the use of our adaptive strategy much harder, since the gradient descent relation ((TJ) relies on the assumption 
that the potential ([6]) is a monotonically increasing function of 7, for 7 moving away from T. However, we have 
implemented our adaptive strategy given by Eqs. (jTJ) , (jT3J) , (JTObJ) , (jTUcJ) and found that it can be successful applied to 
achieve synchronization, provided that 7(0) is not too far away from the stability range T = [1.2,6.5]. Fig. [9] shows 
convergence of j(t) to values in the stability range T for two different choices of the initial conditions, that is, 7(0) = 8 
and 7(0) = 0. 
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FIG. 8: The average synchronization error E obtained from integration of Eqs. JTJ) with 7(£) = 7, versus (constant) 7, for 7 in 
the range [0, 10]. F(x) is the Lorenz equation, h(xi) = x i3 , H = [0, 0, 1], U = 0.75 x 10 3 , t 2 = 10 3 . 
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FIG. 9: We have tested the modified adaptive strategy given by Eqs. © , JT3J , dTObJ , (TUcJ) , for the case of coupled Lorenz 
systems (|14[) , h(xi) = 2^3, TL = [0,0,1]. The sender and the receiver systems are initialized from random points belonging to 
the Lorenz attractor. The adaptive strategy is shown to be successful in two cases, corresponding to two different choices of 
the initial conditions, 7(0) = 8 and 7(0) = 0. Plot (a) shows the time evolution of j(t) from 7(0) = 8, (3 = 10~ 2 , v = 1. Plot 
(b) shows the time evolution of "/(t) from 7(0) = 0, /3 = 10~ 2 , v = 1. 
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IV. CONCLUSIONS 

In this paper we have considered a situation in which a sender system is unidirectionally coupled to a receiver 
and we have studied conditions for the receiver to adaptively evolve the strength of the coupling to achieve stable 
synchronization with the sender. We have proposed a simple adaptive strategy which has been shown to be successful 
in dynamically evolving the coupling strength 7 until it converges to a value in the range of stable synchronization T. 
For cases in which the coupling strength is evolved from initial conditions which are distant from T, we have proposed 
and numerically tested a modified adaptive strategy which includes a saturation on |7|. 

The author is indebted to Prof. Edward Ott for insightful advices and discussions. 

This work was supported by the U.S. Office of Naval Research, contract N00014-07-1-0734. 
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I. INTRODUCTION 

It is known that stability of the synchronized solution of coupled dynamical systems depends on the strength of the 
interaction. Motivated by the observation that real-world networks are characterized by evolving adapting connections, 
some studies have started to take in consideration how the time variability of the strength of the couplings may affect 
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the synchronization dynamics. A few papers have appeared on adaptive strategies, in which the coupling strengths 

(N ; 

, are dynamically adjusted based on information on the states of the systems, e.g., in order to enhance or guarantee 

o; . 

synchronization. In this paper we address the problem of how to evolve the coupling strength from a given initial 

o ■ 

condition to make it converge to a value which corresponds to stability of the synchronized evolution. 



Previous works have addressed adaptive synchronization of chaos. For instance, in [l|, Q], a problem was studied 



in which adaptation was needed in order to achieve synchronization in the presence of external unpredictable events 
affecting the communication between the coupled systems. As an applicative example, consider the case of two sensors 
that seek to synchronize through a continuous signal that they exchange via wireless communication; if an object 
moves across the communication pathway between the sensors, this causes an attenuation on the strength of the 
received signal, which may disrupt synchronization. Therefore, in [l, 2|, adaptive strategies were introduced, having 
the aim of maintaining synchronization with respect to such external (unpredictable) perturbations. 

Here, we will focus on a different but related problem, namely we consider a simple configuration in which a sender 



2 

is connected to a receiver and communication between the two is such that attenuations affecting the received signal 
are negligible; the problem is rather for the receiver to adaptively choose an appropriate coupling form (in what follows 
we will specify our proposed problem to be that of choosing an appropriate coupling strength) with the received signal 
in order to achieve synchronization with the sender. 

To better illustrate this problem, we introduce the equations for the sender and the receiver, 

ii(t) = F(x!(t)), (la) 
x 2 (t) = F(x 2 (t)) + j(t)(H( Xl (t)) - H(x 2 (t))), (lb) 



where x\(t) (^(i)) is the m-dimensional state of the sender (receiver) system, Xi = [xn, Xi2, Xj m ]; F(x) is the 
dynamics of an uncoupled system (hereafter assumed chaotic), F : R m — > R m ; H{x) is an output function, H : R m — > 
R m ; and j(t) is a time varying scalar function measuring the strength of the coupling. Here we assume that H{x) can 
be rewritten as H(x) — Hh(x), where H — [Hi,H.2, ■■■,T~t m ] is a constant m-vector and h : R m — > R is a scalar output 
function. In many practical situations, communication between the connected systems involves only a subset of the 
dynamical state variables of the systems. Therefore, in this paper we consider a situation in which synchronization 
has to be achieved based solely on the scalar signal, h(xi(t)), that system 2 receives from system 1. Our proposed 
problem is to devise an adaptive strategy to dynamically evolve the coupling strength j(t) from an arbitrary initial 
condition 7(0) to a value that is compatible with synchronization. 

As an example of an application for our proposed strategy (to be specified in Sec. II), we consider a sender-repeater 
communication scenario, in which the first device, say a wifi antenna, sends a signal, and the second identical device 
needs to reproduce precisely the same signal to other systems that cannot get the original one. As another example, our 
strategy could be used in applications that take advantage of synchronization of chaos for identification and prediction 
of the dynamics of unknown real systems 3|, U, |5|, |6|. In particular, for the problems addressed in 3|, |4|, |5|, |6j, the 
receiver is a model system which is meant to replicate the dynamics of an unknown true system (the sender) and a 
technique is proposed that aims at identifying the true system parameters by making the receiver synchronize with 
the sender. The technique relies on the choice of a coupling strength which is compatible with synchronization, but for 
this case, an extra difficulty is due to the fact that the dynamical form of the sender's F(x) is unknown at the receiver 
and therefore one is enable to compute the 7-range of stability (i.e., compute the master stability function; for more 
details on this subject the reader is referred to the background section) beforehand. Therefore, if one has to deal with 



such a problem and unless one wants to exclusively rely on a try-and-error approach, it becomes necessary to introduce 
an adaptive strategy to dynamically evolve the coupling strength. We anticipate that the adaptive strategy that we 
present here does not use information from the master stability function and is based exclusively on minimization of 
the squared synchronization error [h(xi(t)) — h(x2(t))] 2 . Therefore, we envision our strategy would be particularly 
useful in applications that use synchronization as a tool to achieve parameter estimation and prediction of unknown 
systems (e.g., for all the applications described in [3j, \4, la, la] ) • 



A. Background 



We note that Eqs. (|la|) and (|lbj) admit a synchronous solution of the type, 



x 2 {t) = Xx{t), 



(2) 



which obeys the sender evolution (JTaJ) . Our proposed problem is to devise a strategy to evolve "f(t) from a given 
initial condition 7(0) in order for the receiver system to synchronize with the sender. Let us now assume that "/(t) 



is constant and equal to 7. Following 7], la llA llll , we know that depending on the choice of the functions F and H, 
synchronization is possible in a certain range of values of 7, say T . To see this, let us linearize the receiver system 
equation 11b[) about @, obtaining, 



5± 2 {t) = [DF(x 1 (t))-*yDH(x 1 (t))]6x 2 (t), 



(3) 



where DF(x\{t)) and DH(x\{t)) are the Jacobians of the functions F and H evaluated about the synchronous solution 
([2]) . Then stability of the synchronous solution depends on the maximum Lyapunov exponent associated with Eq. ([3]) 



(this is a classic result, see 



m 



10( for more details). Once the functions F and H are given, it becomes possible 
to study the dependence of the maximum Lyapunov exponent of (j3]) on 7. The function that associates 7 with the 
maximum Lyapunov exponent of ^ is usually called a master stability function 9|. 



From 



Kj ] and subsequent papers, we know that the range of values of 7, say T, in which the master stability 
function associated with Eq. ([3]) is negative, can be either bounded or unbounded. An example of a bounded range for 
7 is reported in Fig. 1, where for F(x) being the equation of the Rossler system, m = 3, and hixi) = xn, TL = [1, 0, 0], 



the master stability function is found to be negative in the range T = [0.3, 5]. Moreover, we observe from Fig. 1 that 



the master stability function has a minimum at 7 ~ 1.7 (corresponding to the maximum rate of contraction for ([3]) 
and therefore the fastest possible convergence to the synchronized evolution) . Another example of a bounded stability 
range T is reported in [10( for F(x) being the equation of the Lorenz system, m = 3, and h{xi) — x&, TL = [0, 0, 1]. 



2 3 4 5 6 7 

Y 



FIG. 1: Master stability function corresponding to Eq. (3), for our choices of F(x) being the Rossler equation and 
h(xi) — xn, TL — [1, 0, 0]. The zero ordinate line is shown as a guide to the eye (dotted line). 



The problem we address in this paper (a strategy to evolve the s 



■trength of the coupling for 



synchronization) has already been addressed in a number of papers 3 
following adaptive strategies were proposed, 



11 



12 



Q 



achieving stable 



Namely, in 



11 



12, lia], the 



7(t) = a\h( Xl {t)) - h(x 2 {t))\, 
\h( Xl (t)) -h(x 2 (t))\ 



7(*) = 



l + \h(x 1 (t))-h(x 2 (t))\' 



in [12,13], 
in [11], 



(4a) 
(4b) 



where a > is a suitable scalar gain. Equation (|4b[) is approximately equal to Eq. (|4a| for \h(xi(t)) — h{x 2 {t))\ <C 1, 
and is approximately equal to j(t) — a, for \h(x±(t)) — h(x2(t))\ ^> 1. It is worth noting that both strategies are 
based on increasing the coupling strength 7 proportionally with the absolute value of the synchronization error (with 
a saturation in the case of l[4b])). 



We note that the strategies proposed in 



11 



12L 1X31 ] may not be the best suited to deal with coupled dynamical 
systems whose stability is governed by a master stability function as that shown in Fig. 1. For example, for the case 
shown in Fig. 1, if 7(0) is larger than 5, then increasing 7 will not be effective in order to achieve synchronization. 
Moreover, even if 7(0) € T, strategies (|4al4bj) may eventually increase 7 outside the range T (e.g., beyond 5 in Fig. 
1) and yield desynchronization even in a case in which the systems would have synchronized with no adaptation (i.e., 
with 7 = 0). 
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Another adaptive strategy has been proposed in 3[, based on the following adaptation for 7(f), 

7 (t) = -07(f) + - M«a(*))) 2 ]. (5) 

where a > and 5(2;) is a function such that g(x) ~ x for small x and is bounded by a constant C for large x. We 
note from ((SJ, that for large t, 7(f) converges to a _1 <7[(ft.(a;i(f)) — /i(x2(f))) 2 ] and therefore for a monotonely increasing 
g(x), 7(f) increases monotonely with the squared synchronization error. Hence, ([5]) is subject to the same sort of 
limitations as (|4a|4b|) when stability of the synchronized evolution is described by a master stability function as that 
shown in Fig. 1. 

In what follows, we will present a novel adaptive strategy, which is aimed at evolving 7(f) from an arbitrary initial 
condition 7(0) to converge on T. We note that the following limitations may apply to some of the previously reported 
strategies (|4a|4b|5j) . 

1. The strategies (|4al4bl5l) are based on increasing the coupling strength with the synchronization error. Yet, for 
the case of a bounded range of stability T, either increasing or decreasing the coupling strength 7 may be needed 
in order to achieve synchronization. 

2. The strategies (|4al4bl5[) are dependent on the choice of the initial condition 7(0). 

3. The strategies (|4al4bl5[) do not seek to make the coupling 7(f) converge in a neighborhood of its optimal value 
(e.g., the minimum of the master stability function which is observed to be at about 7 = 1.7 in Fig. [T]). 

In what follows, we focus on the common situation that the 7-range of stability is bounded from both below and 
above (e.g., as shown for the case of the Rossler systems in Fig. [TJ. We devise an adaptive strategy that is based 
on both increasing/decreasing 7(f) in order to reach synchronization and will be shown to be independent of the 
initial condition 7(0), as far as 7(0) is not too distant from T. Moreover, though our problem is characterized by a 
continuous of solutions (represented by the range T), our hope is that 7(f) will converge not too far away from the 
value of 7 = 1.7 which corresponds to the minimum of the master stability function, that is to the maximum rate of 
contraction towards the synchronization manifold. 

We wish to emphasize that our strategy described here is only a possible alternative approach to others already 



presented in the literature (e.g., [3 



11 



12 



131 ]) and we assume that many other possible solutions can be found to 
the general problem addressed in this paper. At the same time, we hope our attempt will motivate further studies in 
the same direction, which will provide better and better solutions to our proposed problem. 
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The rest of the paper is organized as follows. In Sec. II, we present our adaptive strategy. In Sec. Ill, we show the 
results of numerical simulations involving our adaptive strategy. In Sec. IV, the conclusions are presented. 



II. ADAPTIVE STRATEGY 

In this section we present our adaptive strategy. We introduce a potential/cost function ^(t), 

*(i)=< [h{ Xl (t)) - h(x 2 (t))} 2 >„, (6) 
where < G(t) > v denotes the sliding exponential average 

j* e -Ht-t')G(t')dt'. From ©, we note that #(t) > 0; 
moreover, due to the chaotic nature of the x's, ^(t) =0 can only be realized if x\{t) and X2{t) are in the synchronization 
manifold, i.e., X\(t) = X2(t). Thus we seek to evolve ^(t) in order to minimize ^(t). To this aim we introduce the 
following gradient descent relation, 

tf . _„« _ „ < lh(xi(m _ >„= WG. (7) 

where (3 > 0, and 

^ = -vG{t) + [h( Xl (t)) - h{x 2 {t))]^^- (8) 

We are interested in how evolves in time. To this aim, we write, 

d dx2 dx2 r dF dH.dx2 rTT/ , xx . , ... , . 

To conclude, our adaptive strategy is fully described by the following set of differential equations, 



7 = 20G, (10a) 
dh T 

G = -vG(t) + [h(x 1 (t))-h(x 2 (t))]— £, (10b) 

0x2 

* = te _ W e + [jff(:Ei W) _ ff(a;2( * ))] - (10c) 



7 



III. NUMERICAL EXPERIMENTS 

In this section we present numerical experiments involving our proposed adaptive strategy (|10p . We integrate the 
set of Eqs. (fTj) . ([TO]) , and we specify our sender-receiver dynamical systems to be described by the Rossler equation, 
in = 3, 



F(xi) = 



— X i2 — Xtf 

xn + 0.2x i2 
0.2 + (x a - I0)x a 



(11) 



1, 2, coupled through h(xi) — xn, TL — [1, 0, 0]. Fig. [2] shows the average synchronization error E defined as, 



E 



(*a -h)pj tl 



\h(x 2 (t)) - h(xi(t))\dt, 



(12) 



from integration of Eqs. Jl} with j(t) = 7, versus 7; the normalization factor p = < (h(xi)— < h(xi) >) 2 > x / 2 , with 
< ... > indicating the time average, is calculated from dynamics of system in the synchronous state, i.e., using 
dynamics from Eq. (|la|) . 




FIG. 2: The average synchronization error E from integration of Eqs. |T} versus (constant) 7, for 7 in the range [—1, 10], 
ti = 0.9 x 10 2 , t 2 = 10 2 , pi = 7.45. 



We first consider the case that 7(0) is not far away from the stability range T = [0.3,5]; namely we consider two 
cases: in the first case we choose 7(0) = 0, and in the second 7(0) = 5.3. In our numerical simulations, system (fl~a|) 
and system (| 1F>|) are evolved from random initial conditions on the Rossler attractor. We choose v — 0.1 so that 
is larger than the characteristic time scale of a chaotic oscillation T ~ 2ir, and (3 to be small enough that 7 changes 
slowly on the time scale of the chaos, that is {3 = 10~ 4 . 
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The results are shown in Fig. [3] (for the case 7(0) = 0) and Fig. [4] (for the case 7(0) = 5.3). In both cases, we see 
that for large enough t, 7(f) converges to T = [0.3,5], e.g., it approaches 2.32 for the case in Fig. [3] and 1.78 for the 
case in Fig. [4] These values are not far away from the minimum of the master stability function at 7 ~ 1.7 (see Fig. 
1). 
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FIG. 3: The upper plot shows the time evolution of Xn(t) (dashed line) compared to X2i(t) (continuous line), while the lower 
plot shows the time evolution of j(t). 7(0) = 0, (3 — 10 -4 , v — 10 -1 . 
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FIG. 4: The upper plot shows the time evolution of xu(t) (dashed line) compared to xn(€) (continuous line), while the lower 
plot shows the time evolution of 7(f). 7(0) = 5.3, /3 = 1CP 4 , v = 10 _1 . 

When 7(0) is distant from the synchronization range T, our proposed problem becomes more complicated. In fact, 
since it may take a long time for j(t) to reach T, X2(t) may move away from x\{t) and the master stability function 
analysis, which describes local stability about the synchronization manifold, may not apply anymore. 

Figure [5] shows a case in which j(t) is evolved from 7(0) = 6.5, for (3 = 1CP 4 and v = 1CP 1 . For this case, 
synchronization is not achieved, and the emergence of an unexpected phenomenon is observed (namely, the emergence 
of a new attractor, corresponding to X2i{t) ~ 0, with both |x23 | and |x22(£)| growing in time in such a way that 
X23(t) — —X22(t)). In particular, from Fig. Owe see that j(t) first decreases from 7(0) = 6.5 to a minimum value of 
about —1.5 and then increases again up to about 44; it is worth noting that 7(i) crosses twice the synchronization 
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range T = [0.3, 5] and neither the first nor the second time it converges to it. We have repeated the experiment in 
Fig. [5] several times, eventually observing a different behavior, that is, j(t) increases indefinitely above 7(0) = 6.5 
and never goes through the synchronization range T (not shown). 

Based on the above observations, we conclude that, in order to improve the effectiveness of our strategy, we need 
to appropriately tune the rate of change of "f(t). In particular, we note that (i) in the case that j(t) is close to T, 
our strategy would benefit from bounding jit) so as to allow enough time for "/(t) to converge in T; (ii) in the case 
that "f(t) is distant from T, our strategy would benefit from choosing large values of (3 so as to quickly adjust ^(t) to 
converge towards V . 
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FIG. 5: The upper plot shows the time evolution of xn(t) (dashed line) compared to X2i(t) (continuous line), while the lower 
plot shows the time evolution of 7(i); in the lower plot, the zero ordinate line is shown as a guide to the eye (dotted line). 
7(0) = 6.5, P = 1CT 4 , v = HT 1 . 



Therefore, we propose to replace -y(t) — 2(3G(t) in (|10a[) by, 



where Eq. (fT"3|) is approximately equal to (|10a[) in the case that \G(t)\ <C 1, while is approximately equal to -y(t) = 
2/3 x sgn(G(t)) in the case that \G(t)\ 3> 1. Note that Eq. (1 13[) corresponds to adding a saturation on 7(f), which 
now is constrained to approximately lye in the range [—13, 0\. 

Figure [S] (Figure [7]) shows the results of numerical simulations in which we have tested our modified adaptive 
strategy, given by Eqs. (fTj). (fT3)) . (|l0b[) , (|10c|) , for a case in which 7(0) = —1 (7(0) = 10). For our experiments in Figs. 
|5]and [71 we have chosen j3 = 10 _1 . As can be seen, in both cases 7(t) is observed to converge to the synchronization 
range T, though it starts from initial conditions that are far away from it. 
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FIG. 6: Modified adaptive strategy. The upper plot shows the time evolution of xu(t) (dashed line) compared to X2i{t) 
(continuous line), while the lower plot shows the time evolution of y(t). 7(0) = — 1, /3 = lO -1 , v = 1. 
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FIG. 7: Modified adaptive strategy. The upper plot shows the time evolution of xii(i) (dashed line) compared to X2i{t) 
(continuous line), while the lower plot shows the time evolution of j(t). 7(0) = 10, (3 = 10 -1 , v = 1. 

A. An experiment with the Lorenz system 

We now consider an example in which our sender-receiver dynamical systems arc described by the Lorenz equation, 
to = 3, 



l0(Xi2 - Xn) 

23x l i - x i2 - XnXiz 



(14) 



n 



i = 1, 2. From [14| we know that when two Lorenz systems (|14l) are coupled as in Jl} with constant coupling j(t) = 7, 
h(xi) = Xi3, TL = [0,0,1], synchronization is stable in a bounded interval of the coupling strength 7, that is for 

7 e r = [1.2,6.5] . 

Figure [8] shows the average synchronization error E, defined in Eq. (|12p. for two Lorenz systems coupled as in 
Eqs. fl}, with constant coupling 7(f) = 7, h(xi) — 2^3, TL = [0,0, 1], as function of the parameter 7; for this case, we 
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17], which 



observe that the transition to the non-synchronous state is characterized by on-off intermittency [15j, 116 
makes the use of our adaptive strategy much harder, since the gradient descent relation ((TJ) relies on the assumption 
that the potential ([6]) is a monotonically increasing function of 7, for 7 moving away from T. However, we have 
implemented our adaptive strategy given by Eqs. (jTJ) , (jT3J) , (JTObJ) , (jTUcJ) and found that it can be successful applied to 
achieve synchronization, provided that 7(0) is not too far away from the stability range T = [1.2,6.5]. Fig. [9] shows 
convergence of j(t) to values in the stability range T for two different choices of the initial conditions, that is, 7(0) = 8 
and 7(0) = 0. 
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FIG. 8: The average synchronization error E obtained from integration of Eqs. JTJ) with 7(£) = 7, versus (constant) 7, for 7 in 
the range [0, 10]. F(x) is the Lorenz equation, h(xi) = x i3 , H = [0, 0, 1], U = 0.75 x 10 3 , t 2 = 10 3 . 
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FIG. 9: We have tested the modified adaptive strategy given by Eqs. © , JT3J , dTObJ , (TUcJ) , for the case of coupled Lorenz 
systems (|14[) , h(xi) = 2^3, TL = [0,0,1]. The sender and the receiver systems are initialized from random points belonging to 
the Lorenz attractor. The adaptive strategy is shown to be successful in two cases, corresponding to two different choices of 
the initial conditions, 7(0) = 8 and 7(0) = 0. Plot (a) shows the time evolution of j(t) from 7(0) = 8, (3 = 10~ 2 , v = 1. Plot 
(b) shows the time evolution of "/(t) from 7(0) = 0, /3 = 10~ 2 , v = 1. 
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IV. CONCLUSIONS 

In this paper we have considered a situation in which a sender system is unidirectionally coupled to a receiver 
and we have studied conditions for the receiver to adaptively evolve the strength of the coupling to achieve stable 
synchronization with the sender. We have proposed a simple adaptive strategy which has been shown to be successful 
in dynamically evolving the coupling strength 7 until it converges to a value in the range of stable synchronization T. 
For cases in which the coupling strength is evolved from initial conditions which are distant from T, we have proposed 
and numerically tested a modified adaptive strategy which includes a saturation on |7|. 

The author is indebted to Prof. Edward Ott for insightful advices and discussions. 
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